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Abstract 

In this contribution we consider sequences of monic polynomials orthogonal with respect to the standard 
Freud-like inner product involving a quartic potential 

{p,q)m= f p{x)q{x)e~^^+'^*^^dx +Mp{0)q{0). 

Jr 

We analyze some properties of these polynomials, such as the ladder operators and the holonomic equation 
that they satisfy and, as an application, we give an electrostatic interpretation of their zero distribution 
in terms of a logarithmic potential interaction under the action of an external field. It is also shown 
that the coefficients of their three term recurrence relation satisfy a nonlinear difference string equation. 
Finally, an equation of motion for their zeros in terms of their dependence on t is given. 

Keywords: Orthogonal polynomials o Freud-like weights o Logarithmic potential o String equation o 
Semi-classical linear functional. 
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1. Introduction 

Let us consider the so called Freud-like inner products 

{p,Q)t= / P{x)q{x)dpt{x), P,(?GP, 

JR 

where dpLt(x) = Wt{x)dx = e~^*'^^'^dx is a positive, nontrivial Borel measure supported in the whole real 
line R. The properties of such sequences of polynomials are known only for certain values of the external 
potential Vt{x). Thus, a well studied example in the literature corresponds to the case Vt{x) = x^ — 2tx^, 
which leads to the inner product 

{P,q)t= f p{x)q{x)e~'^‘'"^'^*'^‘'dx, p,(?GP. (1) 

JR 
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Here t € K+ is a parameter, which, in other contexts, can be interpreted as the “time” variable. We 
denote by {F^}n^Q the corresponding sequence of monic orthogonal polynomials (SMOP, in short) which 
constitutes a family of semi-classical orthogonal polynomials (see (^, [Ig), because Vt(x) is differentiable 
in R (the support of ^t), and the weight function satisfies the distributional equation, known in the 
literature as Pearson equation (see IS) 

[cLii(x)]' = (—4x^ -I- 4te)wt(x). 

where “ ' ” denotes derivative with respect to the variable x. 

We must point out that the class of such a family is s = 2, and that several SMOP related to (HD have 
been studied in the literature. For example, for t = 0 the asymptotic behavior of the corresponding SMOP 
has been deeply analyzed in . In Q a Riemann-Hilbert approach is done to hnd several semi-classical 
asymptotic results for these polynomials. As an application of the semi-classical asymptotics of the ortho¬ 
gonal polynomials, the authors prove the universality of the local distribution o f eig envalues in the matrix 
model with the double-well quartic interaction in the presence of two cuts. In [ig the authors study the 
varying quartic Freud-like weight with V{x) = N -|- for t < 0, where the orthogonality 

takes place on certain contours of the complex plane. They use the Deift/Zhou steepest descent analysis 
of the Riemann-Hilbert problem associated with the corresponding polynomials, to present an alternative 
and a more direct proof of the Fokas, Its and Kitaev result, showing that there exists a critical value for 
t around which the asymptotics of the recurrence coefficients are described in terms of exactly specified 
solutions of the Painleve I equation. In fl9| the author explores the nonlinear difference equation satisfied 
by the coefficients in the three term recurrence relations for polynomials orthogonal with respect to expo¬ 
nential weights, and in the author finds the relation of such Freud’s equations with discrete Painleve 
equations for certain potentials, as V(x) = x^, V{x) = ax'^ + /3x^ or V{x) = x®, among others. The sur¬ 
vey [I^ presents analytic aspects of general orthoranal polynomials associated with several exponential 
weights on hnite and inhnite intervals. The works (9| and [l0| deal with uniform and strong asymptotics 
for polynomials orthogonal with respect to certain exponential weights, and their connections with the 
random matrix theory. Magnus showed in a that the recurrence coefficients of semi-classical orthogonal 
polynomials are shown to be solutions of nonlinear differential equations with respect to a well-chosen 
parameter, according to principles established by D. Chudnovsky and G. Chudnovsky. 

In the linear space P of polynomials with real coefficients, let us introduce the following inner product 


{p,q) = f p{x)q{x)ujtix)dx + Mp{0)q{0) 

Jk 


( 2 ) 


= {p,q)t + Mp{0)q{0), 


with M G R+ . We denote by {Qn}n^o the SMOP with respect to the above inner product. {Qn}n^o is 
said to be the sequence of monic polynomials orthogonal with respect to the measure dv{x) = dp,t{x) + 
M 5q{x), where 5q{x) is the Dirac delta function at x = 0. These polynomials are known as Freud-type or 
Freud-Krall orthogonal polynomials. It is clear that {Qh}n^o is a standard sequence, since the operator 
of multiplication by x is symmetric with respect to such an inner product, i.e. (xp, q) = (p, xg), for every 
p, g G P. 


This inner product was first studied in |l4l | in the particular case t = 0. There, the authors introduce 
a Dirac mass point at x = 0 and provide the holonomic second order linear differential equation with 
varying polynomial coefficients. They also give an electrostatic interpretation for the distribution of zeros 
of the corresponding orthogonal polynomials. 

Following this initial work, the aim of our contribution is twofold. First, we generalize the electrostatic 
interpretation of the zero distribution in terms of a logarithmic potential interaction under the action 
of an external field for values of t greater than zero. Second, we derive an equation of motion for the 
distribution of the zeros of the corresponding Freud-type orthogonal polynomials. 
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The remainder of this Section will be devoted to provide some basic structural properties of {F^}n^o 
to be used in the sequel. 

Proposition 1. Let {F^}„^o denote the sequence of Freud monic polynomials orthogonal with respect to 
m- The following statements hold. 

1. Three term recurrence relation. For every n S N, the recurrence relation for is given by 

(see 0/; 

xF(^{x) = + al{t)Fl_^{x) (3) 

with Ffi{x) = 0 and Fq(x) = 1. An important feature of these polynomials is that the recurrence 
coefficients an(t) in the above three term recurrence relation satisfy the following nonlinear difference 
equation 

^al{t) {al_^{t) + alit) + al+^{t) - f) = n, n ^ 1, (4) 

with ag(<) = 0 and a^(t) = ||cc||^/||l||?. This is known in the literature as the string equation or 
Freud equation (see j /. ]la . (3.2.20)], among others). Moreover, notice that ujt{x) is an even 
weight, i.e. {Fffn'^o is a symmetric sequence of orthogonal polynomials. 

2. We will denote by 

\\FX = {FLKh= [ 

JR 

the corresponding norm. Hence 

(nit) = al{t) ■ ■ ■ al{t) > 0, 

w^th\\Ff\\]=Ut)\\l\\l 

3. Evolution equation. The recurrence coefficients a„(t) in satisfy 

dkit) = akit) (afe+i(t) - al_^it)) , fc > 1, 

where “dot” denotes the derivative with respect to time variable t > 0 (see U (F12)],!M, (24)])- 

4. Structure relation f\T^ . Th. 3.2.1, and p.57j). For every n G N, 

[Fnix)]' = aix,t-n)F(^ix) + b{x,t-,n)Fl_4x), (5) 


where 

aix,t;n) = -Aal{t)x, 

bix, t; n) = ialit) {x^ - t + aUt) + al+^t)) . 

5. Holonomic equation fJTh . Th. 3.2.3, and p.57j). The SMOP {Ff}n^o satisfies the second order 
linear differential equation 


[Fl{x)r + Riix)[Flix)r + Slix)F(,ix) = 0, 


( 6 ) 


where 


Rlix) 

Snix) 


= —4(a;^ — tx') — 


2x 


x^ - t + a],(t) + a(,_^4t) 


= Kit) 


4x^ (^a],_4t) + alit) + aKit) - t - 


- t + alit) + alj^4t) 


-{alit) + aKif) - t) {al_4t) + aUt) - f) + 1 
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The structure of the manuscript is as follows. Section [5] will be devoted to the Freud kernel polyno¬ 
mials. In Section several connection formulas between Freud-type and Freud polynomials are given 
and we briefly analyze the coefficients in the three term recurrence relation for Freud-type polynomials 
and a string-type equation for these coefficients. In Section |31 we provide the lowering and raising opera¬ 
tors associated with Freud-type polynomials, the corresponding holonomic equation, and the electrostatic 
interpretation of the zeros of Freud-type polynomials as positive unit charges interacting according to a 
logarithmic potential under the action of an external field. Notice that the second order linear differential 
equation is obtained as a composition of the ladder operators following the ideas in for some absolutely 
continuous measures and Q for measures with Dirac deltas. Nevertheless, we must point out that in our 
case we consider a measure with an extra time parameter t and a Dirac mass located at a: = 0. The 
explicit expressions of the polynomial coefficients in the holonomic equation are given in our contribution. 
Its role in the study of the electrostatic interpretation is emphasized in this Section. In Section [SJ we 
analyze the zero behavior of Freud and Freud-type polynomials in terms of the parameter t as well as 
we deduce some results concerning the monotonicity and speed of convergence of zeros of Freud-type 
polynomials in terms of the mass M. Finally, in Section [SI we explore some numerical results showing 
the location of certain zeros of Freud-type orthogonal polynomials. 


2. FFeud kernel polynomials 

The n-th degree kernel polynomial 


Kn{x,y,t) = ^ 






(7) 


associated with the polynomial sequence {F^}n^o will play a key role in order to prove some of the basic 
results of the manuscript. Notice that if deg p ^ n, then the n-th kernel polynomial satisfies the so-called 
“reproducing property” 

/ Kn{x,y,t)p{x)wt{x)dx = p{y). 


For X ^ y, according to the Christoffel-Darboux formula, for every n G N we get 


Kn{x,y;t) = 


1 Fi^,ix)Fliy) - Fi^,{y)Flix) 


m\?t 


x-y 

together with the confluent expression 

k^O 

From the above formula 


( 8 ) 


(9) 




Kn{0,0-,t) = 


m\\^ 

[7"n+i]'(0)F^(0) - [FiyiO)F^+M 




and also we conclude that iF„(0, 0;t) > 0, no matter the degree n or parity of the kernel, since it is a 
sum of strictly positive terms. 

The weight function uJt(x) is an even function, so the monic polynomial sequence {Fj^(x)}n^o is 
symmetric (see [ 1 , Ch. I]) and the following expressions hold 


i^^(O) = [F,(]"(0) =0, if n is odd, 

[F;(]'(0) = [F^]'"(0) =0, if n is even. 


( 10 ) 
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From ©, ® and m we have for m > 1 

K2Tn-lix,0;t) = K2m-2{x,0;t) 

1 


K 2 mix, 0 ;t) = 


Fim 


{0)Fi^+Ax) 


( 11 ) 


11^1 


and, therefore, all Kn{x,0;t) has even degree for any n. It is obvious that 

1^2m-l(0, 0; t) = iF2m-2(0, 0; t), 

[FL+iYmLio) 


K2mi0,0;t) = 


( 12 ) 


ll^’l 


2m I \t 


From now on, denote the SMOP with respect to the inner product 


(P:9)t,[fc] = / P{x)q{x) 
Jk 


+2tx 


(13) 


which is a polynomial modification of the measure d^t{x) = e ^ dx called the k-iterated Christoffel 
perturbation. If fc = 1 we have the Christoffel canonical transformation of the measure f^t (see (s^, [3^ 
among others). The SMOP is known as the monic k-iterated Freud kernel polynomials. In 

this particular case, due to the parity of the weight function, when k is odd the perturbed measure is not 
quasi-definite and, as consequence, the corresponding SMOP does not exist. On the contrary, when k 
is even the perturbed measure is still positive definite, and therefore the corresponding fc-iterated Freud 
kernel polynomials are well defined. We will denote by 




k even, 


the corresponding norm. 
t [2l 

Let {Fn (a;)}n>o denote the sequence of 2-iterated monic Freud kernel polynomials, or the SMOP 
associated with the measure dpf^ = x^dpt which is the 2-iterated Christoffel transformation of 
The following two lemmas will be useful in the sequel. 


Lemma 1. The 2-iterated Freud kernel polynomials and the Freud orthogonal polynomials satisfy the 
connection formulas, 


where 


Furthermore, 


where 


X^FY^^Yffx) = i^2m+l(2;) +^2m-l-p2m-l(2^), ^>1, 

Fl'm (a;) = -P' 2 m +2 (a;) + CL ^ 2 m (a;) ) W > 1 , 

,2 _ ll^^2;!fllll?.[2] -[Flm+l]'{^) 

mrn-lWI [FL-l]'{^) ’ 

,2 _ 11 ^ 2 ^ llq[ 2 ] _ -^ 2 ^ 2 ( 0 ) 

WFUM ■ 

a- ^2m-l(a^) = xF2^{x) + [^2m-l ~®2m(0] 1^2m-l(a^)i 
xF2^\x) = ^2*^+1 (a;), 

[^2m-l ~ “ 2771 ( 1 )] > 0 . 


(14) 

(15) 


(16) 

(17) 

(18) 
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Proof. We can expand x‘^Fn^'^\x) in terms of the SMOP as 

n+2 

k=0 

For every fc < n, we have bn^k = 0 by orthogonality. Also we have bn^n +2 = !> since we deal with nionic 
polynomials. The following coefficient 


-'n,n+l 


\X Xn 


WK+iW^t 

vanishes, because the numerator is an integral in the whole K of an odd function, so 


/ CO 

-OO 


Next, we easily see that 


b — £'^ — 

^n,n — Srt — 




[ 2 ] 


m\\^ iiF‘i|2 


t ||2 


> 0 . 


Thus, we have 

x'^F^-^^\x) = F^^^ix) +ilF^{x), (19) 

Shifting the index in the above formula for n odd and n even, we obtain formulas (I14II and (11511 respectively. 
Taking a; = 0 in ([T5)) yields 

^2 _ \\^2’V\\t,[2] _ -Plm+2(0) 

\\FU\^ ~ FUO) ’ 

which is well defined since a; = 0 is not a zero of any Freud polynomial of even degree. In case of Freud 
polynomials of odd degree, the above expression is not defined. Letting a; ^ 0 in (I14p . we may apply 
L’Hopital’s rule obtaining 


F'^ - 

S2m-1 ~ 


l^’2fllll?.[2] _ -[^’2Vl]'(0) 


\\FL-i\?t [J^ 2 ‘,„-i]'(o) ■ 


From dSl) 

and taking a; = 0 we have 


F2m+2ix) = xF^^^^ix) - aL+iW^2m(a;)> 


( 20 ) 


^2m+l 


{t) = 


-^2m+2(0) 

^^2*^(0) ’ 


which essentially equals to ^ 2 m- Combining (fT51) with (Uni we have 

X^F*'J^\x) = XF^^^,{X) + [e2m - 4m+M FL{x) 
= xF^^^^ix). 

Dividing by x we obtain the connection formula (I17F 
Next, from ([31) 

-P2m+i(^) = xF^^ix) - al^{t)F^^_^{x), 
and substituting this expression into ([1^1) gives 

x'^F^'J^l^ix) = xF^^(x) + [^2m-l - «L(0] -P2m-l(a;) 


6 











which yields expression m- The X derivative of m evaluated at x = 0, provides 






( 21 ) 


It is very well known (see [Sj, Ch. I, §8]) that, due to the symmetry of the weight function, there exist 
two n-th degree monic polynomials and such that 

Flm{x) = A^{xy, 

^2m+l(*) = xBm{x ), 

being Bmix) the kernel polynomials, with parameter 0, of An{x). Thus, we have that sign FLiO) = 
signAm{0) = (—1)™ and sign [-F 2 m-i]^( 0 ) = sign Bm-i{0) = (—1)™ so therefore, concerning the sign 
of (I^Tl) we have 

,(0) . A„(0) (-1)^ 


sign 


Fimy 


[FL-ino) (- 1 )’"-^ 

Kim-l ~ ®2m(0] > 0- 


which shows that 


= - 1 , 


This completes the proof. ■ 

t [2l 

Lemma 2. The SMOP satisfies the following three term recurrence relation 

xFf'^'^\x) = Fyfl(x) + alF^yi^x), 

where 


al = 


wFyyiwh 


a^(f) and are given in 


- n-lllt,[ 2 ] 
and Lemma [11 respectively. 


Sn 

£2 ' 
>n—1 


;W- 


t [2l t f2l 

Proof. We can expand xFn’ (x) in terms of {Fn’ }„^o and, for orthogonality reasons, the only terms 
remaining are (x) and (x). The coefficient of Fy'^\ (x) is 1 because we deal with monic polyno¬ 
mials and the other one is 


cti = 


wF^yiwi 


> 0 . 


‘.[ 2 ] 


Let notice that 


Applying (fl^ we get 


Oil = 


{Fy'^‘'\Ff)t,m 


{x^Fy^^\Ff), 


{x^Fl;yiFf_f)y 


= 


{FU2,Fl)t+en{FlFf)t 

{Ff+,,Fi_,)t+en-i{Fl-^,Ff_,), 


Sn M-*- nil 


i ||2 
t 


en-i\\Fi.y\i 

a 


^2 

>n—1 


\(t). 
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This completes the proof. ■ 

Let Xn,k = Xn,kit), fc = 1,...,n, be the zeros of F^{x), Fn^‘^\x), respectively, arranged 

in an increasing order. All of them are real and simple. By parity reasons, these zeros are symmetrically 

[21 [2] [2] [2] 

arranged with respect to the origin. That is, —xl^\ = Xn'n, 2 = „_i and so on. 

i [ 2 I 

We next prove that the zeros of G 2 m{x) = xF 2 ^_-^{x) and F^^lx) interlace. Concerning the zeros 
g 2 m,k = k = 1,... ,2m, of G 2 mix) are the same zeros of F 2 ^\-^{x) except one more zero at 

the origin, so G 2 m{x) has a double zero at a: = 0. Thus, we have g 2 m,i = ^ 2 m-i u ^ ~ 1; ■ • ■ with 

[2] [2l 

®2m-l,m = 52m,m = 52m.m+l = 0 and g2m.r = a^2m-l.r-l; r = m + 1, . . . , 2m. 

Both G 2 m{x) and F^^^x) are even polynomial functions, so their respective graphs are symmetric 

with respect to the point x = 0, which means that we only need to prove interlacing in the positive real 

semi-axis, being the situation in R_ the reflection with respect to the y-sods. Hence, we only need to prove 

t [ 2 I 

that, for X > 0, between two consecutive zeros {x 2 m,k , X 2 m,k+i) there are only one zero of F 2 ^_i{x). 
Consider m evaluated at the positive zeros of Fl^{x) (i.e., X 2 m,r > 0, with r = m -I- 1,..., 2m). We 
have 

X2m,r F 2 ^_l(x 2 m,r) = [^2m-l ~®2m(0] ^ 2 m—lix 2 m,r), 

and taking into account (fTSl) . we get 


sign[F 2 ’J^}_^{x 2 m,r)] = sign[F 2 ^_-^^ix 2 m,r)], r = m -k 1,..., 2m. 


( 22 ) 


Thus, from iSTIh . the relation between the zeros of G 2 m{x) and T' 2 )|fli(x), the symmetric reflection with 
respect to the y-axis, and the well known fact that the zeros of T 2 m-i(^) interlace with the zeros of 
^ 2 m{x), we obtain the following interlacing property 

Theorem 3. The inequalities 

X2m,l ^ 52m,1 ^ ^2m,2 ^ ' * * *\ X2m,m ^ 52m,m = 0 = ^2m,m+l ^ ^ X2m,2m—1 ^ 52m,2m ^ ^2m,2m ; 

hold for every m € N. 


3. Connection formulas 

We can expand Qh{x) in terms of the SMOP {F^}n^Q, which is an orthogonal basis of P, as follows 


n —1 


Q*n(x) = F^ (x) 4- Xn,^F*{x), 


2=0 


where 


. _ {Fl (x) ,Qlix))^ 

'^n.i — 


't nt , , , ^ 

n,2 — I I I |9 7 0 ^ i ^ Tl 1. 


\\Fm 


From m and d?]) the above equality becomes 

Qlix)= F*^ (x) - MQi (0)iL„_i(x, 0; t), 
and evaluating the above expression at x = 0, we deduce 

Fm 


(23) 


Qm = 


l + MiL„_i(0,0;t)' 


(24) 




Hence, from ([3]) and (H]) , (1331) reads as 

xQnix) = K+lix) + A^F^ix) + 


where 


At 


B 


t 

n 


II 

Ir ri-1 

al{t) + ■ 


MF^,{0)F^,{0) 

MjFmf 


/ l + MKr,{0,0;t) \ 
{l + MKn-ii0,0;t)) ■ 


Remark 1. Due to the fact that uJt{x) is an even weight function, x = 0 is always a zero of Ff(x) for n 
odd. Then, Ff(0)Ff_i{0) = 0 for every 1, and therefore A^ = 0 for all positive integer n. 


From the above remark, we have 


xQtix) = F^^^ix) + BlF^_^{x), (25) 

Q2n+l(x) = ^2n+l(^)' 


Introducing the notation 


we get 


1 + Mi4:„(0,0;t) 
l + MKn-i(0,0;t) 

Bn = al(t) bl. 


(26) 


This yields an expression for the ratio of the energy of polynomials Q\i(x) and Ff(x) with respect to 
the norms associated with their corresponding inner products. 


Proposition 2. Let H-H^ be the squared norm of Freud-type monic polynomials with respect to (0). Then 

\m\^ ^ l + MKn{ 0 , 0 ;t) ^ 

||F*||2 l + MKn-i{0,0;t) - ■ 

Moreover, = 1 when n is odd and > 1 when n is even, i.e., for every to ^ 0, 


HQ 


2m+l I 


\\Fl 


2m+l 


||2 
lit • 


Proof. Taking in account 

llQnf = {Qi{x).X-) 

= {Qii.x),Fl(x)) 

= {Qi{x),Fl(x))^ + MQi(Q)F^^{Q) 

and using (l24ll . we get 


IIQ: 


t ||2 


t ||2 




M{F^nW 


t ||2 
t 


imi 

\\F% 


1 + MKn-i (0,0; t) 

(l + MjF„_i(0,0;t)) + 


1 + Mif„_i(0,0;t) 


tii2 1 + Af(0,0; t) 

l + MiF„_i(0,0;t)' 
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Evaluating m at a: = 0 we have 6^ = 1 for n odd and 6^ > 1 for n even. This gives the result when 
combined with the above equation. ■ 

As a consequence, we get 


Theorem 4. Let be the sequence of monic Freud-type polynomials orthogonal with respect to 

(0]. Then 

Qlm+iix) = Fi^+i{x), m > 0, 

MF^O) 


Q2mix) = F^^{x) - 


1 + MK2m-i{0, 0; t) 


K 2 Tn-iix, 0 ;t), mfil. 


Next, we provide an alternative way to represent the Freud-type polynomials of even degree Q 2 m{x) 
in terms of the polynomials F^^^x) and the 2-iterated monic Freud kernel polynomials Fn^‘^\x). This 
representation will allow us to obtain the results about monotonicity and asymptotic behavior (presented 
below in this work) for the zeros of Q 2 mix) in terms of the parameter M present in ([2]). We only need 
to consider Freud type polynomials of even degree, because they are the only ones affected by variations 
of M. 


Theorem 5 (Connection formula). The sequence {Q 2 m}m>o represented as 


Qlmix) = F 2 ^{x) -I- MK2m-i{0,0;t)G2m{x), m > 1, 


(27) 


with G 2 Tnix) = xF 2 J^}_j^{x), Qlrnix) = H 2 mQnix) and K2m = 1 + MK 2 m-l{ 0 , O', t) > 0. 


Proof. Let K 2 m be given by the positive quantity K 2 m = 1-1- MK2m-i{0,0',t). From the expression for 
Qimix) in Theorem m we have 

K2TnQlmix) = K 2 mFl^{x) - M Fl^(Q)K 2 m-l{x ,Q]t) , m > 1. 

Being K 2 m > 0, we call Q 2 m{x) = K 2 mQ 2 m{x) the polynomial with the same zeros that Q^mix)- Next, 
dl) and (fT^ yields 

Qlmix) = F^^{x)M [K2m-ii0,0;t)F^^{x) - F^^{0)K2m-iix,0;t)] 


= FL(x)-M 


{FL-iy{o)FLio) 


\\Fl 


2m— 1 \ \t 


FLix) - 


FLjO) Fi^_,{x) 
[FL-ino) X 


— F2^{x) + MK2Tn-li0,0',t) 


FLix) - 


^2m(0) 


[J^2‘^-l]'(0) 

On the other hand, replacing m into dm) and dividing by x yields 

A.[2] ^2m(0) ^2m-l(^) 


xFff^Liix) = F2 ^{x) - 


[^’2‘™-l]'(0) O: 


(28) 


as 


From the above expression, we can finally rewrite 

Q\m{x) = Fl^{x) + MiL2™-i(0, 0; t)xFl'j^\,{x). 


This completes the proof. ■ 

In the remaining of this section we will focus our attention on the coefficients of the three term 
recurrence relation satisfied by the Freud-type SMOP. Since Q*n{x) are standard and symmetric, they 
satisfy the following fundamental recurrence relation. 
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Proposition 3. The polynomials Qn(x) satisfy the three term recurrence relation 


xQiix) = Qi+iix) + jri{t)Qi-i{x), (29) 

where 

7„(t) = ^al{t). (30) 

Proof. We expand xQ\{x) in terms of the SMOP {(5(j}„^o3'iid, taking into account (HJ, the Lemma 
follows. ■ 


Proposition 4 (String equation). The coefficients i30\) of the above three term recurrence relation for 
{Qn}n>o so.t'isfy the following nonlinear difference string equation 


47^(0 


bi-2 

bi 




bi 


7n(i) + ^r^7n+i(^) - ^ j =n. 


Proof. It is enough to replace (l30ll in the string equation ([4]), and then the Proposition follows. 


4. Holonomic equation and electrostatic model 

Next, we give details of the second order linear differential equation satisfied by {Qh}n^o when t > 0. 
First, from ([3]) we can rewrite (|2^ as 

= Ai{x,t;n)Ff{x) + Bi(t;n)Ff_j^{x), (31) 


with 


Ai{x,t',n) = X, 

Bi(t;n) = a^(t) (bl, - 1) . 

In order to obtain the ladder operators and the second order linear differential equation, we follow a 
different approach as in (isl . Ch. 3]. Our technique is based on the connection formula (I23|) . the three 
term recurrence relation ([3]) satisfied by the SMOP {Ff}n^o, and its corresponding structure relation 

®. 

We begin by proving several lemmas which are needed for the proof of Theorem [9l 
Lemma 6. For the SMOP {Qh}n^o emd {Ff}n^Q we have 

AQhi.^)]' = Ciix, t;n)F^{x) + Di{x,t;n)F^_^{x), (32) 

where 

Ci{xA-,n) = -Aal{f) [hl^x^ + (6(, - 1) {al_^{t) + al{f) -t)] , 

Di{x,t;n) = 4al(t) X {al+j^(t) + bl^ [x^ - t + al{t)]) - ^al{t){bl^ - l). 

The coefficients Ai(a:,t;n), Bi(t;n) are given in h31\) . b*^ is given in F2()\) . and a{x,t]Ti), b{x,t',n-) come 
from the structure relation W- 

Proof. Shifting the index in ([5]) as n —>• n — 1, and using ([3]) we obtain 

[Fn-iix)y = a{x,t]n)F^{x)+ b{x,t;n)F^_-i^{x), (33) 
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where 


a{x,t;n) = -4:[x'^ - t + al_^{t) + al{t)] , 

b(x,t;n) = Ax [x'^ — t + . 

Next, taking derivatives with respect to the variable x in both sides of m. we get 

x‘^[Q\{x)]' = x‘^[F^{x)]' +xBi(t;n)[F^_-i^{x)]' - Bi{t;n)F^_j^{x). 

Substituting ([5]), (1551) and the expression for Ai(x,t;n), Bi{t;n) into the above expression the Lemma 
follows. ■ 

From Proposition!^ = 1 when n is odd, so the above results can be simplified as 
x[Q\m.+i{x)\' = C'i(a:,t;2m + l)F 2 *^+i(a;) + Di{x,t]‘2.m + l)Fl^{x), 

where 

C7i{a;,t;2TO +1) = -4a|^+i(t) 

Di {x, t-2m+l)= (t) x (x^ - t + (t) + (t)^ , 

which corresponds to (O. 

Lemma 7. The sequences of monic polynomials {Qh}n^o o,nd {Ff}n^Q are also related by 

xQ*n-i{x) = A2{t-,n)Fl{x) + B2{x,t;n)F^_^{x), (34) 

x[Ql^_-i^{x)]' = C2{x,t;n)F^{x) + D2{x,t;n)F^_;i^{x), (35) 

where 

A2{t]n) = 1 
B2{x,t;n) = 

and 

C 2 {x, t] n) = (&(,_! - 1) (Axalit) + j) - 4x6(,_i [x^ - t + + al{t)] , 

D2{x,f,n) = Aalif) (l - [x^ - t + al{t) + + Ax%l^_^ [x^ -t + ai{f)\ . 

Proof. The proof of (1551) and (1551) is a straightforward consequence of (|5]), ([5]), (15T|) . and Lemma [S] ■ 
By Proposition [2] it is obvious that (154)) is exactly the first equation of Theorem |4] if n is even and 
(1551) can be simplified as 

x[Q\m-i{x)]' = C2{x,t]‘2m)F^^{x) + D2{x,t]‘2m)F^^_j^{x), (36) 

where 

C2{x,t]2m) = -Ax{x'^ - t + ajrn-lit) + > 

L> 2 (a;,t; 2m) = 4a;^ (a;^ - t + a 2 „(t)) . 

An equivalent formulation of (|56l) is (O if we substitute FL (x) in the above relation according to ([5]) 
when the index n is shifted by 2m — 1. 

The following lemma shows the converse of (15T1) - (155)) for the polynomials Ff{x) and Ff_^{x). Indeed, 
we express these two consecutive polynomials of the SMOP {Ff}n^o in terms of only two consecutive 
Freud-type orthogonal polynomials of the SMOP {Qh}n^o- 
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Lemma 8. For t > 0, 


xhi_iF^{x) = xhl_-^Q\^{x)-al{t){b\,-l)Qi_-^{x), 

xbl_^Fl_^{x) = (6^_i - 1) Ql{x) + xQl_j^(x). 


(37) 

(38) 


Proof. Note that (PTl) and (IM)) can be interpreted as a system of two linear equations with two polynomial 
unknowns, namely F^(x) and F^_-^^{x). Hence from Cramer’s rule, we have 

A(x, <; n) = xHl_^ + al{t) {b*^ - l) {bl^_^ - l) . 

From Proposition we know that b*^ = 1 when n is odd and b*„ > 1 when n is even. Therefore, the 
product {b\ — 1) {b\_i — l) is always zero, and the Lemma easily follows. ■ 

Theorem 9 (Ladder operators). Let a„ and be the differential operators 

d 


Cln — 


a), = 


A 2 (t;n) Ci{x,t;n) 
B 2 {x,t]n) Di(x,t',n) 

di(x,t;n) C 2 (x,t]'n) 
Bi{t;n) D 2 {x,t',n) 


+ A(x, t] n) 


dx' 


- A{x,f,n) — , 
dx 


satisfying 


^n[Qni^)] — 
^i[Qh-i{x)] = 


Ai{x,t','n) Ci{x,t]n) 
Bi{t;n) Di{x,t]n) 

A2{t;n) C2{x,t\n) 

B2{x,t',n) D2{x,t',n) 



(39) 

Qiix). 

(40) 


Let point out that all the above expressions are given only in terms of the coefficients in (EB, ©, 
m, and (1351) . Thus, a„ and ajj are, respectively, lowering and raising operators associated with 
the Freud-type SMOP {Qhfn^o- 

and (1551) one 


Proof. The proof of Theorem [5] follows from Lemmas [5H51 Replacing (I571) - lia»ll m 
obtains the ladder equations 


^ , _ Ci{x,t;n)B 2 {x,t;n) — A 2 {t;n)Di{x,t;n) 


= 


A(x, t\n) 


Q„{x) 


Ai(x, t; n)Di{x, t] n) — Bi(t; n)Ci(x, t; n) . 


and 


( , C2{x,t-,n)B2{x,t;n) — A2{t;n)D2{x,t;n) ^ 




A{x, t] n) 


A{x, t; n) 

Q*nix) 


Qn-lix) 


Ai{x, t; n)D 2 {x, t; n) — Bi{t; n)C 2 {x, t; n) . 


A(x, t] n) 

which are equivalent to (l39l) and (l40l) . This completes the proof of Theorem |9l 


Qn-l{x), 


Theorem 10 (Holonomic equation). The Freud-type polynomial Q*„(x) satisfies the holonomic egua- 
tion (second order linear differential eguation) 


A{x,t;n)[Ql{x)]” + B{x,t;n)[Ql{x)]' +C{x,t;n)Ql,{x) = 0, 


(41) 
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where 


A(x, t] n) 
B{x, t; n) 

C(x, t; n) 


A{x,t;n) [W A(a;,i;n)} + i; n) [^' 2 ,i(x, n) - 


A{x,t;n)W 4'2,i(a;, <; n)} + t; n) 




'i'i, 2 ix,t\n)]) , 
^'i,2(a;,<;n) 
^'2,2(a;,<;n) ’ 


with 


'i>ij{x,t;ri) 


Ai(x,t',n) Cj{x,t',n) 
Bi(x,t',n) Dj(x,t',n) ’ 


i,j = 1,2. 


Moreover, the polynomials A{x,t;n) and C{x,t;n) are even functions and the polynomial B{x,t;n) is an 
odd function, whose coefficients are showed in Table [7J 


Proof. The proof of Theorem [TU] comes directly from the ladder operators given in Theorem [HI The 
usual technique consists in applying the raising operator to both sides of the equation satisfied by the 
lowering operator, i.e. 


4 [Qi-iix)] = al 


1 


Thus, 


^i,2(a^,t;n) 
^ip{x,t;n) 
which becomes 


aniQni^)] - A{x,t;n) 


^ip{x,t;n) 
d 


^n[Q ni^)] 


dx \dti i(x, t;n) 


= '^2,2ix,t;n)Ql,{x). 
an[Q^,(x)]) ='i’2,2ix,t;n)Ql,{x) 


'}’i^2{x,t-,n)'Sr2,l{x,t-,n) / XI ^'i,2(a:,t;n)A(a:,t;n) 

-A{x,t-,n)£ [ l"fl%'Xn) [Qu{x)\ + xi/t.uS) IQU^)]') = 'f2.2{x,t;n)Qi{x). 


The intermediate computations above yield 


d f'^2,i{x,t]n) . 


dx \t;rr) 


mx)] = 


± (^Mi^LiqUx)]'] = 


4'i_i(x,<;n) ’l>' 2 ;^(x,t;n) 

'I>i,i(x,t; n) A(x,t;n) 
tl/[j^(x,t;n) A'(x,t;n) 


Q^(x) 


+ 


'I' 2 ,i(x, t;n) 




[Qn{x)]', 


mx)y 


A{x,t-,n) ^ 


and 


4'i,2(a;,t;n)4'2,i(a;,t;n) 


^'ni(a;,<;n) 


[Q\i{x)\ + 


^ 1 . 2 ( 2 :, ^;?^)A(a::,^;r^) . 


[QUx)]", 


'^iA{x,t;n) 


A(x, t; n) 
'i’i,iix,t;n) 
A{x, t; n) 


^i^i{x,t;n) 5'2,i(x,t;n) 

4''2_i(a::, t; n) 

t; n) A{x,t;n) 
'^[ilx,t;n) A'{x,t;n) 


-mx)]' 

A{x, t, n)'^t 2 ,i{x, t; n) 


A'^{x, t; n) 


Unix)]' - 


Unix)]'' 




isix,t',n) 

= '^2.,2ix,t-,n)Qy{x). 

Combining all the above expressions, and after some cumbersome computations. Theorem 1101 follows. 


Corollary 1. For every nonnegative integer n, when n is odd (1411) is eguivalent to Aix,t-,n) times ([HI) 
and, in the other case the polynomial coefficients of (1411) are contained in Table 
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Proof. The procedure is to observe that in the case where the degree is odd, m is reduced to the 
second order linear differential equation ([6|), which is satisfied by T 2 m+i( 2 ^) since 

B{x, t; 2m + 1) = A[x, t; 2m + 1) R^m+ii^) C{x,t', 2m + 1) = A{x, t; 2m + 1) Sl^_^_i{x). 

On the other hand, if the degree is even then we get & 2 m-i = 1- ■ 

For a deeper discussion of ladder operators we refer the reader to [isL Ch. 3]. We next provide the 
second order linear differential equation satisfied by the SMOP {Qn}n^o taking into account the measure 
/it is semi-classical. This is the main tool for the further electrostatic interpretation of zeros. Once we 
have the second order linear differential equation satisfied by the SMOP {Qhjn^o it is easy to obtain an 
electrostatic model for their zeros. We will study the asymptotic behavior of the position of the movable 
constant charges involved in the external field. As we have shown in section [3] Freud-type orthogonal 
polynomials and Freud polynomials of odd degree coincide. Thus, in this Section we shall derive the 
electrostatic model for the zeros in the case when n is even. 

Let us evaluate (IdTIl at the zeros {y 2 m.i(t)}?™i of the polynomial yielding 

[^2m (2/2m,2 (t))] ,(t),t;2m) 

[Q 2 m {y2mAt))]' A(y 2 m.»(t) , t; 2m) 

+ 8y2m,»(t) blm " t) + alm+lA)] 

4^m2/2m.iW +4//L,zW[(a 2m(t)^m ^2771 Aim 

-^ + 4(yL.i(t) -1 y2m,i{t)) > 

y2m,7{t) 

where 

Am = (4aL(t) [almit) + Am-lit) - t] “ l] “ l) {Am “ l) ■ 

The above equation reads as the electrostatic equilibrium condition for {j/ 2 m,i(t)}i™i. Having 

u{x,t;2m) = +4a;2([aL(t)&L -t]&L + aL-Hi(t)) + (42) 

the previous condition can be rewritten as 


2m - 

y _ i _ 

- y 2 m,iit) 


_l_ 1 M'( 2 / 2 m,i(t),t; 2 m) 
2 u[y 2 m,iit)A'A'm) 


1 

y2m,i{t) 


+ 2(?/L,i(t) -ty2m,i(t)) =0, 


which means that the set of zeros {j/ 2 m,i(t)}i™i are the critical points (zeros of the gradient) of the 
total energy. Hence, the electrostatic interpretation of the distribution of zeros means that we have an 
equilibrium position under the presence of an external potential 

V^^\x) = ilnu(a;,t;2m) 

where the first term represents a short range potential corresponding to a unit charge located at the real 
zeros of the quartic polynomial u{x,t',‘2m), and the second one is a long range potential associated with 
the Freud weight function. 


5. Zeros of Freud-type SMOP 

Having in mind the techniques shown in (l6| , we will first study the motion of zeros of time depending 
of the polynomial F^x) (Theorem [TT]) and, finally, in subsection 15.21 we give the differential equation 
that the zeros of the polynomial QAx) satisfy. 
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5.1. Equations of motion for zeros of {Ff{x)}n^Q 

Notice that the three term recurrence relation (|31) implies that 

+ [al^^{t) + al{t)\Fl{x) + al{t) al_^{t) Fl_.^{x) 
dF^ 

for all n ^ 1 and Ff^^x) = 0. The requirement on —^{x) is that 


(43) 


QF* 

n 

dt 


dt 

n—2 


(x) = 


i=0 


with \>n,i = 


££n pt 

dt ’ 


\\m: 


■. From 


/ OO 

F^{x) Fl{x)ujt{x)dx, O^i^n — 2, 

-OO 


differentiating the above expression with respect to t, we obtain 


0 = 


apt 8F* ' 

F!{x)^{x)+F^^{x)^{x) 


uJt{x)dx + 2 / x^ F^{x) Fl{x)ujt{x)dx. 


dFf 


Hence, since (x) is a linear combination of the first n elements of the sequence {Ffcjfe^Oi using (H51) 
we get 


\ dt ' 


Fi ) — 2a„(t) a„_]^(t) ||-F'„_ 2 |L Sn- 2 ,i, 0 ^ i ^ n 2, 


and 


dF* 

dt 


(x) = -2alit) al_:^_{t) F^_ 2 {x) 


(44) 


as claimed for t > 0 and n ^ 1 with Ffi{x) = 0. 
We can now formulate this result 


Theorem 11 (Equation of motion for zeros of Ff{x)). Let n be a positive integer and t > 0. If 

Xn,i{t),... ,Xn,n{t) are the n zeros of F^(x), then 


dxn,k{t) 

dt 


i.k (t) 


(^1/cW - ^ + flnW + a^+1 W) 


Proof. Given n ^ 1 and t > 0, we have 

Fl^{xn,k[t)) =^, l^kifn. 

If we differentiate the above equation with respect to time t, we get 

dxn,k{t) dFf 


dt dx 


fx) 


^=Xn,k{t) 

On the other hand, from (O and (I45L we obtain 
dF* 


dF* 

+ ^(x„,fe(t)) =0. 


(45) 


(46) 


dx 


(x) 


^ — ^n,k (^) 


(^) ; ^5 1 (O) * 
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Hence, using (HU) . (HSl) and the above formula, it follows that 


dxn,k{t) 

dt 


K-2{Xn,k{t)) 


2 i^lkit) -t + al{t)+ Fl_^ {xn,k{t)) 


Finally, we obtain the result since the three term recurrence relation ([3]) yields 

K-2{xn,k{t)) = Fl^_^{Xn,k{t)) ■ 


( 47 ) 


5.2. Equations of motion for zeros of Freud-type orthogonal polynomials 
It is required that 


dt 


iy) = Qliy), 




with bn ,2 = 
relations 


oQn 
dt ’ 


IIQ 


,tl|2 


. We first compute the coefficients b„,i, 0 ^ ^ n — 2. From the orthogonality 


0 = 


Qn{y)Q\{y)^t{y)dy + MQI^{Q)Q\[0), O < f < n - 2 . 


Taking the partial derivative with respect to the variable t in the above equation we obtain, for 0 ^ z ^ 

n — 1, 

'^{y),Q\{y)) = -2{QUy),y^Ql{y)),- 

Applying ([8]) and (f23|) in the above equality we get 

bn,z — d; 3, 


Since 


(Qliiy),y^Qiiy))t = (Puiv),y^Ql{y))t - mqI{ 0) {Kr,-i{x, O; t), y^Qliy))^. 


Then, we get 


where 


dQl 

dt 


iy) — ^n,n—2 Qn— 2 ( 2 /); 


bn,n—2 — 


bl,-2 




M[Ff{0)]^ 


(1 + Mif„_i(0,0;t)) \\Ff_^\\J 

Let us denote by yn,i{t), ■ ■ ■ ,yn,nit) the n zeros of polynomial Ql^{y). From (1^51) it follows that 

r yn,kit)bi-2 ( u^\ 

(yn,fc(i)j — bn,n-2 ^2 fj.-^Qn-l\yn,k{t)) ■ 

We can now state the analogue of Theorem [TT] for the zeros of Freud-type polynomials. 


(48) 
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Theorem 12 (Equation of motion for zeros of Ql^{y)). Let n be a positive integer and t > 0. If 
yn.iit ),... are the n zeros ofQl^{y), then 


where 


(Iyn,k 

dt 


(t) = 


~yn,ki^) '^n,n-2-yr—^-^JT—^Jp^Qn-l{yn,kit)') 
Ck,i{n,t)Ff{yn,k(t)) + Ck, 2 {n,t)Ff_-^(yn^k{t)) ’ 


Ck,iin,t) = l-4:al{t)[bl^ylj,{t) + {bi-l){al{t)+al_i{t)-t)], 

Ck,2{n,t) = ialit) ynMt) ^ W) +a^+i(i)] • 

Proof. Given n ^ 1 and t > 0, we have 


Qh{yn,k{.t)) = 0, l^k^n. 

If we differentiate the above equation with respect to t, we get 


dyn,k 

dt 


it) 



y^yn,k{t) 


dQi 

dt 


{yn,k{t)) = 0- 


(49) 


(50) 


We only need to compute the y derivative of Qn(y) and to combine (H51) with the above expression in 
order to get the result. Differentiating (1^ and applying ([3]), (O, and ((431) . we have 


= Ck,i{n,t)Fl{yn,k{t)) + Ck,2{'^,t)F^_^{yn,k{t)) ■ 

V=Vn,k{d 

Hence, combining (H51) and the above formula with (1501) the Theorem follows. ■ 

In the case when n is odd, the above Theorem provides (|47ll . which establishes the formula of Theorem 

El 


/,N dQ\ , . 
yu,k{t) -^{y) 

dy 


5 . 3 . Behavior and monotonicity with M of the zeros o/Q|„(a:) 

Let assume that yn,k, k = 1,2, ...,n, are the zeros of Qh{x) arranged in an increasing order. From 
the analysis done before, it is clear that the zeros 2 /„ ^ when n is odd are not affected by the mass M. 
Next, we analyze the behavior of zeros y2m,s = y2m,s{M), s = 1,..., 2m, as a function of the mass M 
and we obtain such a behavior when the positive real number M goes from zero to infinity. In order to 
do that, we use a technique developed in [j, Lemma 1] and (ill . Lemmas 1 and 2] concerning the behavior 
and the asymptotics of the zeros of linear combinations of two n-th degree polynomials /ira, 5 n S P with 
interlacing zeros, such that f{x) = hn{x) + cgn{x), c ^ 0. From now on, we will refer to this technique as 
the Interlacing Lemma. Here the linear combination of two polynomials of the same degree 2m is given 
by ([27)) . and ^ 2 m play the role of hn{x), gn{x) respectively. 

In order to apply this technique, we need to show that the hypotheses of the Interlacing Lemma are 
fulfilled. First, in Theorem[3]the interlacing of the zeros of and G2m was proved. In our computations, 
we will only deal with the zero behavior in the positive real semi-axis, because the behavior in K_ 
follows by reflection through the y-axis by symmetry reasons as usual. Thus, from (1371) . the positivity of 
^2m-i(0, 0; <), and Theorem|3]we are in the hypothesis of the Interlacing Lemma, and we immediately 
conclude the following results about monotonicity, asymptotics, and speed of convergence for the zeros 
of Q2m,i^) ill terms of the mass M. 

Let us define the monic polynomials 

GI^(x) = x{x - g2m.l)ix - g2m.2) ■ ■ ■ (x - g2Tn,m-l) 

= x{x - 41 - 1 , i)(a; - 4I-1.2) ■■■{x- 41 - 1 ,m-l) 
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and 


^m(^) — x(^X (J 27 n,m+ 2 ){x 92m, 771 + 3 ') {x 9‘2m,2m') 

= x{x — X 2 m-i^m+l)i^ ~ ^ 2 m-l,m+ 2 ) ' ' ' {x — a^ 2 m-l, 2 m-l)) 
such that C^^ix) = Gl,{x)G^^{x). 

Theorem 13. In the negative real semiaxis, the following interlacing property holds 

3^2?71,1 y2m,l 92m,1 ^ ^2m,2 y2m,2 <•••*< 9‘2m,m—l ^2m,m ^ y2m,m ^ 92m,m — 0- 

Moreover, each = yn,i{M) is an increasing function of M and, for each I = 1,... ,m, 

lim y 2 m,l{M) = g 2 m,l , 

M—¥ 0 C) 


well 


1- Tirr ] _ '^2m(52m,/) 

lim M[g2m,l — y2m,l\ — T 7 /n rt-T-Wm T' 

xV2m —1 (O 5 0? [^mj (y2m,z) 


Applying symmetry properties through y-axis, in the positive real semiaxis the following interlacing pro¬ 
perty holds 

0 — g2m,771+1 ^ y2m,m+l ^ X2m,m+1 ^ g2m,m+2 ^ * * * ^ X2m,2m—1 ^ g2m,2m ^ y2m,2m ^ X2m,2m- 

Moreover, each yn,r = yn,r{AI) « decreasing function of M and, for each r = m + 1,..., 2m, 

lim y 2 m,r — 92771,r 7 
M —¥00 


as well as 


lim M[y 277 i,r - 92 m,r] 
Af —>00 


^277li9277l,r) 


K277i-i{0, 0; t)[G^]'{g 2 m,r) ' 


Notice that the mass point at a; = 0 attracts two zeros of Q 277 i{x), i.e. when M -+ oo, it captures 
y2m,m and j/ 2 m,m+i at the same time. 


6. Numerical experiments 

We next provide some numerical experiments using Mathematica® software, dealing with the zeros of 
Freud-type polynomials More specifically, we will show the position of the two symmetric and 

closest-to-the-origin zeros of some even polynomials of the sub-sequence • We choose Q\{x) for 

the following first experiment varying t. We show the location of its zeros for several values of t and M, 
and we also show the position of the source-charges of the short range potential 

Vshortix) = ^\nu{x,t;2m), m = 1,2,3,..., 

which are the zeros of the polynomial u{x, t; 4) defined in (H^ . In tables[3im[5]and|n]we provide numerical 
evidence of the position of its zeros when t is equal to 1/2,1, 3/2, and 2 respectively, for several choices 
of M. Notice that the polynomial u{x, t; 4) has exactly degree four, and its zeros are always two real 
and two simple conjugate complex numbers. We also remark that we recover the results in |l4l . Cor. 3.5] 
when t = 0. 

Figure [T] illustrates the change in the even Freud-type polynomials when M varies as described in 
Theorem m We enclose the graphs of Q\{x) for three different values of M. The black continuous, 
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dashed, and dotted lines correspond to M = 0, M = 0.2, and M = 0.6, respectively. We also include, 
with different tones of gray color, the graphs of Q^x) and Ql{x), showing that the odd degree polynomials 
are not affected by the variation of the mass M. 

Finally, the last two tables [3 and [5] show the position of the zeros of Freud-type polynomials of Ql{x) 
and Q\q{x) and the zeros of the corresponding ghost polynomials. Notice that the zeros of the ghost 
polynomials continue to be two real and two complex conjugate numbers. 
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Figure 1: The graphs of Q\(x) and Q\(x) (gray) and Q\{x) for some values of M (black lines). 
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7 
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Table 1: Coefficients of the polynomials A(x^ t; n) and B(x^ t; n) for every integer n. 
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Coefficient 
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4ai„(t)(oi„+i(t) + - t]) 

A{x, t; 2n) 
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2 

0 
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11 

-16oi„(t) 


9 
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7 

-4ai„(t){l - 4tai„+i(t) + 4o^„(t)[5*„ - 1]^ + + 4t26‘^ 
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-4tai„(t)(l -4oi„(t) [a^nit) + - t] [bln “ l] ) (^2n “ 0 


3 

-2ai„(i)(l - 4ai„(t) [ai„(<) + oi„_i(t) - t] [&|„ - l]) (6|„ - l) 


1 
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Table 2: Coefficients of the polynomials A(x, t; 2n) and B(x^ t; 2n) for every integer n. 


t = 0.5 



M = 0 

M = 0.002 

M = 0.05 

M = 0.5 

M = 5 

M = 10 

M = 50 

< 34 ( 2 ;) 

± 1.1640 

±1.1639 

±1.1623 

±1.1516 

±1.1318 

±1.1286 

±1.1257 


± 0.4839 

±0.4836 

±0.4755 

±0.4154 

±0.2256 

±0.1689 

±0.0794 

u{x, t; 4) 

±0.94861i 

±0.94869* 

±0.9505* 

±0.9618* 

±0.9827* 

±0.9863* 

±0.9898* 


0 

± 0.0144 

±0.0689 

±0.1528 

±0.1219 

±0.0942 

±0.0455 



Table 3: Zeros of Q 4 '^(x) and 'u(x,0.5;4) for some 

values of M. 






t = 1 





M = 0 

M = 0.002 

M = 0.05 

M = 0.5 

M = 5 

M =10 

M = 50 

Qlix) 

±1.3002 

±1.3001 

±1.2988 

±1.2891 

±1.2659 

±1.2615 

±1.2570 


±0.6156 

±0.6153 

±0.6084 

±0.5533 

±0.3335 

±0.2551 

±0.1227 

u{x, t; 4) 

±0.7653i 

±0.7654* 

±0.7687* 

±0.7894* 

±0.8354* 

±0.8454* 

±0.8563* 


0 

±0.0171 

±0.0827 

±0.1962 

±0.1800 

±0.1424 

±0.0703 


Table 4: Zeros of Q\{x) and u{x, 1; 4) for some values of M. 
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t = 1.5 

M = 0 

M = 0.002 

M = 0.05 

M = 0.5 

M = 5 

M = 10 

M = 50 

Ql(x) 







±1.4485 

±1.4484 

±1.4474 

±1.4395 

±1.4120 

±1.4047 

±1.3964 

±0.8059 

±0.8057 

±0.8010 

±0.7603 

±0.5363 

±0.4290 

±0.2176 

u{x, t; 4) 







±0.5175i 

±0.51791 

±0.52631 

±0.57141 

±0.6631 

±0.68991 

±0.72821 

0 

±0.0233 

±0.1125 

±0.2770 

±0.3021 

±0.2477 

±0.1261 

Table 5: Zeros of Q\^{x) and u(x, 1.5; 4) for some 

values of M. 





t = 2 




M = 0 

M = 0.002 

M = 0.05 

M = 0.5 

M = 5 

M = 10 

M = 50 

Qiix) 

±1.60437 

±1.6043 

±1.6038 

±1.5989 

±1.5717 

±1.5594 

±1.5406 

±1.0429 

±1.0428 

±1.0408 

±1.0220 

±0.8736 

±0.7644 

±0.4490 

u{x, t; 4) 

0 

±0.06051 

±0.19061 

±0.33161 

±0.45641 

±0.48851 

±0.57481 

±0.1487 

±0.1613 

±0.2539 

±0.4256 

±0.5625 

±0.5117 

±0.2832 


Table 6: Zeros of Q^{x) and u{x, 2; 4) for some values of M. 






t = 1 





M = 0 

M = 0.002 

M = 0.05 

M = 0.5 

M = 5 

M = 10 

M = 50 

Qlix) 


±1.51614 

±1.51612 

±1.5153 

±1.5103 

±1.5018 

±1.5005 

±1.4993 


±1.0730 

±1.0729 

±1.0711 

±1.0600 

±1.0403 

±1.0374 

±1.0346 


±0.4530 

±0.4526 

± 0.4445 

±0.3846 

±0.2044 

±0.1524 

±0.0714 

u{x^ t; 6) 


±0.91641 

±0.91651 

±0.91851 

±0.93001 

±0.95011 

±0.95331 

±0.95641 


±0 

±0.0139 

±0.0665 

±0.1444 

±0.1109 

±0.0852 

±0.0409 


Table 7: Zeros of Qq{x) and u{x, 1; 6) for some values of M. 
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Qioi^) 


±1.79469 

±1.79467 

±1.7942 

±1.7921 

±1.7896 

±1.7893 

±1.49286 

±1.49284 

±1.4922 

±1.4888 

±1.4849 

±1.4845 

±1.17419 

±1.17414 

±1.1730 

±1.1674 

±1.1608 

±1.1600 

±0.7931 

±0.7930 

±0.7907 

±0.7789 

±0.7647 

±0.7630 

±0.2950 

±0.2947 

± 0.2858 

±0.2291 

±0.1067 

±0.0780 

±1.1107* 

±1.1108* 

±1.1117* 

±1.1166* 

±1.1223* 

±1.1230* 

±0 

±0.0110 

±0.0509 

±0.0950 

±0.0589 

±0.0440 


Table 8: Zeros of Q\q(x) and u{x, 1; 10) for some values of M. 


±1.7890 

±1.4841 

±1.1593 

±0.7616 

±0.0359 

±1.1236* 

±0.0206 
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